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Abstract 

We study the 2D Ising model on three different types of lattices 
that are topologically equivalent to spheres. The geometrical shapes 
are reminiscent of the surface of a pillow, a 3D cube and a sphere, 
respectively. Systems of volumes ranging up to 0(10^) sites are simu- 
lated and finite size scaling is analyzed. The partition function zeros 
and the values of various cumulants at their respective peak positions 
are determined and they agree with the scaling behavior expected 
from universality with the Onsager solution on the torus (i^ = 1). For 
the pseudocritical values of the coupling we find significant anomalies 
indicating a shift exponent 7^ 1 for sphere-like lattice topology. 
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1 Motivation and introduction 



Most of the Monte Carlo studies of spin models have been done on hypercubic 
lattices with periodic boundary conditions, i.e. with the topology of a torus. 
Recently, however, there has been a growing interest in the effect of other 
topologies. Some of that interest is motivated by the advance of quantum 
gravity models, some by problems related to the role of topological excitations 
of the particular system. 

Various investigations [1-6] have been carried out on lattices topologically 
equivalent to the surface of a hypersphere in D+1 dimensions. On sphere-like 
surfaces a loop can be continuously contracted to a point: The fundamental 
homotopy group is trivial. This has immediate consequences for the dynamics 
of string-like objects like monopole loops (in 4D pure gauge theory with 
U{1) gauge symmetry) or the boundaries of clusters in 2D spin models. We 
therefore might expect a different approach to the thermodynamic limit and 
different finite size corrections to scaling. 

At first order phase transitions the phase mixture obtained in the ther- 
modynamic limit may be influenced by local changes, like flxing one spin, or 
changing the boundary conditions. The critical exponents however, should 
not depend on such effects. At second order phase transitions one expects 
that the critical properties are not affected by boundary conditions or by 
whether one works on lattices with torus or spherical surface geometries 
if either becomes flat in the inflnite volume limit. Such an assumption of 
universality should be tested in practical examples, which is one of the mo- 
tivations for this work. Another question is whether the flnite size scaling 
(FSS) ansatz is general enough to persist. The important leading terms and 
the corrections to them may have different size depending on the lattices and 
their topology as has been observed in the 4D and 2D studies. 

The Ising model is explicitly solved for torus geometry [7j even on flnite 
lattices |Hlini- It has a well understood 2^*^ order phase transition. FSS of 
the bulk quantities is dominated by the leading term and provides a good 
example of the power of this method to determine critical indices. 

Here we want to present a Monte Carlo study for the Ising model on 
different 2D lattices with the topology of the surface of a sphere: the surface 
of a 3D cube, a pillow-like structure and the cubic surface projected to a 
sphere. In all these cases we determine cumulants and partition function 
zeros on different lattice sizes in order to study the size dependence and FSS 
in detail. After the introduction of these geometries in sectEl the simulation 
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and multi-histogram analysis are discussed in sect El and the results in sectEJ 
We conclude that there are sizable differences to the usual results on torus- 
like lattices, but that universality prevails. Preliminary results of this work 
have been presented in [H]. 

2 Lattice geometry 

The most straightforward approach to constructing a sphere-like lattice would 
be triangulation. We want to stay as close as possible to the original action, 
however. This is partly motivated by the related studies P and 0, where 
one wants to keep a 4-link plaquette structure. Also, two of the lattices 
studied here can be viewed at as combinations of the usual square lattices 
glued together at the boundaries; thus in the thermodynamic limit even the 
non-universal critical coupling ought to agree with that for torus geometry. 

The torus is our reference geometry; in addition three other lattice ge- 
ometries were used for our simulations (cf. fig. Q). 

Torus TO [N]: The solution for the Ising model on an x lattice with 
periodic b.c. is known jHl E] and we will use this as the standard to 
compare our results with. 

Pillow PI [N]: This is the surface of a cube with N x N x 2 sites, where A^ 
will be called the base length in the subsequent discussion. Basically 
it is made out of two N x N lattices glued together at the edges. The 
curvature is concentrated in the 4x2 corners, elsewhere the lattice is 
locally fiat. This lattice PI [N] has 2N'^ sites and 4(A^2 - 1) links. 

(Dual) Cube SH [N]: This lattice is dual to the surface of a 3D cubic lattice 
of base length A^. Each plaquette of the cubic surface is identified with 
a site of SH . This has the advantage, that each site has 4 links to 
nearest neighbors, like for the torus. The lattice has 6(A^ ~ 1)^ sites 
and 12(A^ — 1)^ links. The curvature is concentrated on the 8x3 corner 
points. 

Sphere S [N] : That lattice is defined by the sites on the surface of a A^^ 
cube, projected onto the unit sphere. The site-link connectivity struc- 
ture is identical to that of the cube's surface, the only difference lies 
in weight factors in the action, to be discussed below. The lattice has 
6(A^ -1)2 + 2 sites and 12(A^ - 1)^ links. 
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The total energy (or action, if used in the context of quantum field theo- 
ries) for the model is given by 

E = -J2'Wx,flSxSx+fj, (1) 

X,fl 

where the sum runs over all links (at sites x in directions fi). The spin 
variables Sx have values G Z(2). The weight factors Wx^p. are equal to 1 
for the torus and the lattice types PI and SH . The partition function is 
determined by the sum 

Z = J2exp{-(3E{C)) (2) 
c 

over all spin configurations C. 

For lattices PI and SH the curvature is concentrated around the corners. 
The deficit angle denotes the deviation of the sums of angles of plaquettes 
(or triangles) at a given site from the flat-space value 27r. For PI its value 
is 7r/2 on each of the 4x2 corners, for SH it is ir/G on each of the 8x3 
corners. It vanishes on all other sites, as it does on the torus: The lattices 
are fiat almost everywhere. The total curvature on the sphere-like lattices 
therefore has the value Air or Euler number 2, as compared to for the torus. 
Euler's relation for sphere-like lattices is Ug + np — ni = 2 (n^, Up, and ni 
denote the total number of sites, plaquettes and links, respectively) whereas 
for torus topology this sum is zero. 

Both lattice types PI and SH can be imagined as built out of fiat N x N 
pieces, glued together along their boundaries. The deviation from the torus 
shaped lattice thus disappears at least as fast as a boundary contribution 
0{1/N). The contribution from the corners, where the curvature is concen- 
trated, is suppressed 0{1/N^). We expect that, although not a universal 
quantity, the value of the critical coupling jSc (in the thermodynamic limit) 
coincides with that for the torus. 

In order to mimic a truly spherical lattice more closely, we also took the 
structure of a cubic surface lattice projected to the sphere as shown in fig. [T] 
Since the link connectivity structure does not change, the difference has to 
be expressed by modifying the weight factors accordingly. Christ et al. ITU] 
discuss a possible definition for triangulated (random) surfaces, which obeys 
certain properties, that are necessary for the consistency of the continuum 
field theory at a second order phase transition. 
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Our lattices are built out of quadrangles and therefore we have to modify 
that method for the derivation of the weight factors. We start with the scalar 
continuum field theory with the Euclidean invariant kinetic term 

Jd'x{d,<f){x)f . (3) 

(In fact, we can define the Ising model as a discretization of a scalar field 
theory in the limit of infinite quartic self coupling.) The discretization on 
a lattice replaces the derivative term by the nearest neighbor differences, 
d^(j){x) — >• {(px+ji — (px)/lx,x+fi, where lx,y denotes the distance between the 
neighboring two points. The area differential d'^x has to be replaced by an 
area element A^^fi assigned to each lattice link (or site term). The discretized 
kinetic term then leads to the link term of the lattice action (^, with the 
weight factor Wx,ii = 2^3,,^//^^^+^. 

The total lattice volume has to be split into lattice area elements at- 
tributed to the link contributions. These lattice area elements are constructed 
with help of the dual lattice, where each site may be defined as the barycenter 
of the plaquettes of the original lattice (cf. fig. EI). Links of the dual lattice 
are drawn between sites corresponding to neighboring plaquettes. Each link 
on the original lattice corresponds exactly to one link of the dual lattice. 

For simplicity we make two approximations: 

• We assume that all quadrangles are flat; in reality the four corners are 
not co-planar. In the thermodynamic limit the lattice becomes locally 
flat and the error vanishes. 

• We neglect angular distortions and approximate each quadrangle by 
a rectangle; this allows us to write the area assigned to a lattice link 
X as the product of the lengths of the dual link with that of the link 
itself, Ax j:^ = ^ I'^j^ x l^ f^. The error introduced due to this approx- 
imation does not vanish in the thermodynamic limit and introduces 
slight distortions from the regular spherical surface. However, the cur- 
vature is still smeared out over the whole lattice, although somewhat 
non-uniformly. As will turn out later, that appears to be no serious 
problem for the finite size analysis. 

^In |10| the center of periphery was used; there is no such point in a general quadrangle, 
therefore the center of mass seems to be a reasonable alternative. 
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This way the total lattice volume is distributed over all links. 

With these approximations we find that the weight factors for the lattice 
geometry ^[A^] are just the ratio of dual link length to the link length 

w^^,A = ^x,ii/kf^ (4) 

equivalent to the factor suggested in [THj for triangulations. 

For this particular sphere-like lattice we do not expect, that the thermo- 
dynamic limit value of the critical coupling is identical to that of the torus. 
Since we have introduced a change of the action affecting all links, (3c will be 
renormalized. Still we expect universal behavior of the critical exponents, to 
be checked numerically. 

The computer programs use index tables to deal with the lattice geometry. 
For the lattices 5'[A^] the weight factors are precalculated and tabulated. In 
our discussions we will refer to the lattice volume 

V = T.^x,i. (5) 

X,fl 

as the typical size quantity. For lattices TO , PI and SH this is just the 
number of links. A length scale may be defined as L = \JV/2, such that its 
value is for a torus TO [N]. 

In another study of the Ising model [3j regular honeycomb lattices folded 
to tetrahedron shapes have been used. There it was possible to obtain series 
expansions as well as Monte Carlo results. In [l] the Ising model was studied 
with Monte Carlo methods on triangulated random lattices with sphere-like 
topology. We compare these results with ours in the section HI 



3 Simulation details 

In the Monte Carlo simulation we studied lattices with the base length A^ = 
16, 32, 64 and 128. (The torus results were determined up to A^ = 256.) The 
lattices have the volumes given in tabled 

For the updating of the spins we used the Swendsen-Wang cluster al- 
gorithm pi]. For each lattice size the bulk energy histograms were deter- 
mined for various values of the coupling and then combined with help of the 
Ferrenberg-Swendsen multihistogram technique This leads to an opti- 
mal estimator for the distribution densities Pl{E) for the partition function 

ZL{P) = Y.PLiE)expi-pE) (6) 

E 
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Base Length 


TO 


PI 


SH 


S 


16 


512 


1020 


2700 


2724.8 


32 


2048 


4092 


11532 


11650.0 


64 


8192 


16380 


47628 


48147.0 


128 


32768 


65532 


193548 


195726.9 



Table 1: The volumes V as defined in the text; the equivalent linear extent 
L = {V/2)^ ranges from 16 (for the torus) up to 313 (for the sphere). 



and allows to determine various moments {f{E)). Of course the statistical 
accuracy deteriorates when one evaluates these quantities outside the domain 
of /3-values, where one determined the histograms. Suitable overlap between 
the individual histograms is necessary for a reliable application of the method. 
We determined only the energy histograms and analyzed only observables 
in the even sector of the model and not the observables that involve the 
magnetization. 

We measured the specific heat, the Challa-Landau-Binder cumulant JH] 
and another 4*^^ order cumulant suggested by Binder (cf. the review [T^): 

cviP^L) = ^{{E-{E)r), (7) 
VclbW^i-) - ^^^p , (,«) 

The positions and values of their respective extrema are used for the FSS 
analysis. 

Eq. (jH)) defines implicitly an analytic continuation to complex values of 
f3 not too far away from the real axis. Therefore it is possible to determine 
the nearby zeros of the partition function jT3] in the complex /3-plane, the 
so-called Fisher zeros jTHj. As will be demonstrated below, in particular the 
imaginary part of the zero closest to the real axis provides a high quality 
estimator for the critical exponent u with small corrections to the leading 
FSS behavior (cf. ^7] for a recent high statistics study of the Ising model in 
4D, where it was possible to identify the logarithmic corrections to scaling 
on basis of the Lee- Yang singularities ^Hl)- So real and imaginary parts of 
the closest Fisher zeros provide further (even) observables. 
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For each lattice size we simulated the system at up to 20 different values 
of (3 between 0.41 and 0.47. The integrated autocorrelation times for the 
energy varied between r ~ 3 — 7. For each size we produced between 2 
and 5 x 10^ independent configurations. The errors were estimated with the 
jackknife algorithm, i.e. from the variation of the results for the analysis of 
subsamples of the raw data. 

Throughout the discussion of the data we use L = \/V/2. The fit quality 
is expressed through the goodness of fit parameter Q. □ 



4 Results and analysis 

4.1 Finite size scaling and topology 

From the usual scaling hypothesis [IHl HHI one expects for the singular part 
of the free energy density the scaling behavior 

/(r,L) = L-i/^/(rLV^l), (10) 

where r = (1 — P/Pc) denotes the reduced coupling and L is the length scale. 
From this one derives the scaling behavior of the cumulants. At a second 
order phase transition we expect (for D = 2) 

C (L) ~ l^'^^^ (for«>0) 

^max[-L^) - j oQnL) (fora = 0) ^ ^ 

T/ fr\ { L'^'"-^ (fora>0) 

VcLB,mUL) - I ^-,1^^ (fora = 0) ^^2) 

/ 0(1) + 0(L-/-) (for«>0) 
Ui,minKi^} - j 0(1) + 0(1/^^)2 (fora = 0) ^ ' 

\mzo{L) ~ L^^'" (14) 
l3,{L)-(3, ~ L-^ (15) 

with Josephson's law a = 2 — Dv. The asymptotic value of t/4 depends 
on the details of the distribution density p{E) and is e.g. 3 for a Gaussian 
distribution. For the Ising model the Onsager solution gives u = 1 and a = 0. 



^The goodness of fit parameter Q is defined as Q = r(^^^^, x^)/r(^^^^), where n is the 
number of fit points and p is the number of fit parameters. It is the integrated probabihty 
over all larger than the measured one. 
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We denote by (3c{L) our definitions for pseudocritical points: Tlie posi- 
tions of tlie extrema in tfie cumulants. The so-called sliift-exponent A is for 
many models equal to 1/z/, but not necessarily so in general; this relation is 
not a necessary conclusion of FSS (cf. the discussion in [IHl)- We return to 
this issue in the discussion of the results. 

A priori we know nothing about the absolute size of the multiplicative 
coefficients in the scaling formulas. They depend on the details of the lattice 
geometry and topology and on the boundary conditions. 

The FSS behavior comes from the rescaling properties of the bulk quan- 
tities. The effect of changing the boundary properties may be responsible for 
further contributions. The non-homogeneous distribution of the curvature 
in our lattices PI and SH might also be responsible for additional (constant) 
terms in the total free energy. 

Actually, since the total curvature is an invariant, there may be another 
contribution, which — relative to the bulk contribution OiV) — becomes 
irrelevant in the thermodynamic limit. It has been shown 20^, that the total 
free energy for finite 2D systems with non-singular metric and smooth bound- 
aries has at criticality (in addition to boundary terms 0{L)) an asymptotic 
contribution proportional to InL. The proportionality constant is a prod- 
uct of the central conformal charge and the Euler number (vanishing for the 
torus). Thus this contribution depends only on the topology of the system, 
not on the shape of the boundary. 

All these contributions to the free energy are suppressed 0(1/L^) or 
0(lnL/L^) relative to the leading term. In the absence of a strict theory 
we therefore might expect corresponding additive corrections terms in the 
FSS relations (fTT|) - (ffe|). 

4.2 Cumulant values and partition function zeros 

The values of the cumulants at their respective pseudocritical points provide 
information on the critical exponents according to (fTTj) - (fT^ . As discussed, 
they may have geometry dependent corrections. However, in our data we 
find qualitatively excellent agreement with the scaling of the torus-results 
and no significant indication of geometry-corrections. 

Fig. Elshows that the specific heat scales with InL, as expected for a = 0. 
A power law fit gives a value for a compatible with and has a larger x^'- 
The logarithmic behavior is preferred. 

Comparing the results for the higher order cumulants Vclb and f/4 we 



9 



also find excellent agreement with the torus results, if compared at the corre- 
sponding scales L, and with the expected scaling behavior for a = according 
to (fT^ - (fT^ . In particular the results for Vclb he on top of a common curve 
for all geometries. 

The finite size dependence of the positions of the partition function zeros 
confirms this observation. In fig. |3] the nearest Fisher zeros for sphere- 
like lattices are compared with those for toroidal lattice. The real part is 
substantially closer to the thermodynamic limit. Its scaling properties are 
discussed below together with the pseudocritical points derived from the 
cumulants. 

The imaginary part of the closest Fisher zero appears to profit from the 
smallness of the deviation of the real part from the thermodynamic value. 
The log-log plot (fig. demonstrates the excellent scaling signal and a fit 
of the form 

Im ZoiL) = aL-^/" (16) 

gives u = 0.9964(46) (goodness of fit Q = 0.39) for PI , u = 0.9975(57) 
(Q = 0.77) for SH and u = 1.0023(54) (Q = 0.12) for S lattices. For all 
these lattices the result is in perfect agreement with the value z/ = 1 of the 
toroidal lattice. 

We conclude, that the imaginary part of the first partition function zero is 
an optimal observable for extracting the critical exponent u. It appears to be 
least affected by correction to scaling due to lattice topology and boundary 
effects. 

In this light the excellent scaling behavior of the specific heat should not 
be too surprising, since the peak value is directly related to the vicinity of 
the closest Fisher zero. Since the partition function is proportional to the 
product of all zeros, 

Z^lli/3- Zi) (17) 

i 

the specific heat includes for its singular part the contribution 

The closest zero therefore contributes a term oc (Im zq)^'^ to the peak value 

of Cy 

The fact, that the closest zero approaches the real axis (with increasing 
lattice size) almost perpendicular is also clearly exhibited by the shape of 
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the specific heat itself. In fig. IHlwe compare the torus results with those for 
the cubic surface lattices for equivalent lattice volumes. The approach to the 
infinite volume case (Onsager solution) is in a much more symmetric way 
than for the torus lattice. 

4.3 Pseudocritical points 

We discuss here the pseudocritical values derived from 

- the peak positions of the specific heat, 

- the minima positions of the other cumulants Vclb and f/4, 

- the real part of the position of the closest zero in the complex /?-plane. 

For the lattice geometries PI and SH we expect (see the discussion in sect. 
12)) that in the thermodynamic limit the critical values jSc coincide with those 
of the torus, and we therefore present these results in direct comparison. 
For the lattice type S the asymptotic value of the critical coupling will be 
somewhat different and we discuss these results separately. 

4.3.1 PI and SH lattices 

It turns out, that both lattice geometries have very similar behavior and agree 
(except for the smallest lattice PI [16]) even numerically with each other, if 
compared at corresponding volumes. 

As fig. [7| clearly exhibits, there is an obvious difference in the FSS behav- 
ior compared to the usual torus results. The overall size of the corrections 
to the thermodynamic value of the critical coupling are much smaller for the 
sphere-like lattices. The leading FSS behavior of f3c{L) for large L should 
follow (fT3j) . For the Ising model on a torus the shift exponent isA = l/z^ = l 
jnj. This leading behavior linear in 1/L is evident in the figure. However, for 
PI and SH another effect seems to blur this picture: A possible (but clearly 
very small) linear term is dominated by contributions nonlinear in 1/L. 

As discussed in the leading linear term may vanish even for x M 
torus geometry, depending on the ratio N/M. In particular it vanishes in the 
limit M — * 00, where the leading behavior becomes O (In L/L^) jU]. There are 
also other specific models and cases, where \^ l/v ^Hj- As mentioned below 
(I15|l also the topology may give rise to additional terms [20] in the free energy 
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(per unit volume) proportional to the Euler number and to 0(lnL/L^); it is 
unclear how these affect the pseudocritical points in our particular situation. 

This observation, that the dominating behavior appears to be non-linear 
in 1/L, was also made in a study of the Ising model for a honeycomb lattice 
on a tetrahedron surface Both, series and Monte Carlo results led to a 
value A = 1.745(15) ISj- The value of v obtained there from the correlation 
length and the specific heat agreed with the Onsager value. 

We therefore fit our data for the pseudocritical points to 

(3,{L)-(3, = aL-' + bL-'^ , (19) 
pc{L)-l3c = aL^^ + bL-HnL , (20) 

with the Onsager value for /3c = | In (1 + V2). 

The data for each of the 4 definitions of pseudocritical coupling (from 
cy, VcLB, and Re Zq) appears to be consistent for both geometries PI 
and SH . We therefore use one set of parameters a and b different for each 
definition but identical for the two geometries. The value of A is assumed 
to be universal for all definitions and both geometries. We fit to data for 
N > 32. The results according to ()19j) are given in table El and are plotted 
in fig. [7| they fit the data perfectly with a x^/d-f. ^ 1.2 (Q = 0.26). 



Lattice 


A 


Par. 


cv 


VcLB 




Re Zq 


PI,SH 


1.76(7) 


a 


0.002(4) 


0.000(13) 


-0.011(5) 


0.009(10) 






b 


-0.83(18) 


-4.29(96) 


-0.37(13) 


-0.45(24) 


S 


1.71(10) 


a 


0.028(10) 


0.025(28) 


0.032(11) 


0.029(11) 






b 


-0.62(13) 


-3.34(86) 


-0.24(14) 


-0.29(13) 



Table 2: Scaling law coefficients a, b and exponent A for the various lattice 
geometries. Fits are according to ()19|1 to data for = 32 . . . 128 for PI and 
SH and = 16 . . . 128 for S 

As expected from looking at the data we find only a small contribution 
to the term 0(1/L), compatible with zero for all observable except for f/4. 
Removing this term altogether seems conceivable, although is quadrupled 
in this case. The second term clearly dominates. However, the resulting 
value A = 1.76(7), although consistent with the results for the tetrahedron 

is not stringent. In fact, allowing for the second ansatz (EIHl give almost 
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the same fit quality and would be indistinguishable in the figure. Also fixing 
A to a value 2 is still compatible with the data. 

We conclude that we are in a situation where a possible leading FSS 
term 0(1/ L) has an (almost or completely) vanishing coefficient and the 
subleading terms dominate. This resembles the Ising model on a cylinder 
with infinite extension in one direction. It cannot be decided, whether the 
corresponding term is of form (fTn|) or ()2()|1 . 

4.3.2 Spherical surface lattices 

Now we turn to the approximate spherical surface topology (S). Whereas 
for the pillow and cubic surface lattices the curvature is concentrated in 8 or 
24 points here it is more or less uniformly distributed among all sites of the 
lattice. The total curvature (the Euler number) remains constant. 

In fig. ISlthe peak positions of the specific heat and the other cumulants 
are plotted together with fitted curves according to (fT!?|) . Here Pc is a free 
parameter, otherwise we follow the procedure discussed above, i.e. one com- 
mon value of A but different parameters a and b depending on the observable. 
Since we have fewer data we include the data from the smaller lattices with 

= 16. This is also justified by the overall smaller deviations from the 
asymptotic value. 

The resulting values are also given in table |21 We find a behavior in 
agreement with the other sphere-like lattices. The contribution of the term 
0(1/ L) is again very small and the non-linear term dominates again. The 
overall variation of the pseudocritical points with L is for most observables 
smaller than for the other lattice geometries. The fitted critical temperature 
ispc = 0.43883(3) with a xV^^-/- - 0-5 {Q = 0.81). The value of A = 1.71(10) 
is consistent. 

Again the fit to ()20|1 gives results of comparable quality and corresponding 
curves would be indistinguishable in the figure. 

5 Conclusion 

We have performed a high statistics Monte Carlo study of the Ising model 
on lattices of various size and different shapes, all with sphere-like topology. 
Our FSS analysis led to the following conclusions. 



13 



We have found explicitly, that the Ising model on a spherical surface 
topology lies in the same universality class as the planar Ising model 
with periodic boundary conditions — topologically the surface of a 
torus. Our results demonstrate that universality holds, independent of 
the lattice geometry. This agrees with similar conclusions obtained for 
tetrahedral lattices [Sj and random lattices |lj of sphere-like topology. 

However, some observables are not well suited to find the expected 
leading FSS behavior. Different observables vary in their sensitivity. 
Of the studied quantities (in the even sector of the Ising model) we find 
that the imaginary part of the Fisher zero closest to the real axis has 
the smallest (in fact: not identifiable within our accuracy) deviations 
from the leading FSS behavior. Related to this quantity, also the peak 
value of the specific heat scales according to the FSS formulas with the 
Onsager values for the critical exponents without further (identifiable) 
corrections. The values of the other cumulant have larger statistical 
errors but are also in agreement with the torus results. 

The change in the topology class infiuences the size of the FSS contri- 
butions. This appears to affect in particular the pseudocritical points. 
We find no significant contribution of 0(1/L), which is the dominant 
FSS term for the torus pseudocritical points (with a shift exponent 
X = l/u = 1). Instead we find that the FSS behavior is dominated by 
a term 0{L^^) with a mean value A ~ 1.74(6) (averaging the results for 
PI , SH and S). A compatible value was obtained in an independent 
study for tetrahedral lattices [3]. The behavior seems to be universal 
for all sphere-like lattices, independent of the details of the geometry. 
This contribution is also consistent with a term 0(lnL/L^); such a 
term describes the FSS of the specific heat peak position for the Ising 
model in cylinder geometry [H]. It also has been argued, that a term 
of that kind contributes to the free energy per unit volume for systems 
with non-zero Euler number [20]. Unfortunately this implies that none 
of these observables is qualified to derive the critical exponent u. 

In general we find that the studied sphere-like lattices have smaller 
corrections to the infinite volume behavior than one observes for the 
torus (i.e. periodic boundary conditions). The approach to the ther- 
modynamic shapes is faster and in a more symmetric way. 
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• Our analysis of the Fisher zeros of the partition function is consistent 
with this picture. With increasing size the closest zero approaches the 
real axis almost perpendicular. (Note, that this behavior is a finite size 
behavior and is not identical to the asymptotic impact angle - defined 
e.g. as the angle between the first and second zero). The results for the 
spherical surface geometry appear to be closest to the thermodynamic 
behavior in general. For different lattice (sphere-like) geometries the 
FSS behavior is consistent if one chooses the size variable L = 
where V is the number of links. This choice appears to be preferable 
over the base length A^. 

In general our conclusion are consistent with other results on sphere- 
like lattices for the Ising model jHl 1^ and with similar observations in other 
models [H], also in higher dimensional lattices [2]. 

How can one explain the more symmetric and seemingly faster approach 
to the thermodynamic limit, that occurs for the sphere-like lattices as com- 
pared to the torus? One may argue, that on the 2D torus there are two glob- 
ally distinguished directions. As soon as the correlation length ^ approaches 
some fraction of the linear size, the system notices the loss of its rotational 
invariance. Larger clusters may then span in the two distinguished directions 
of the lattice. On sphere-like surfaces (although there is local orientation) 
there are no globally distinguished directions. Rotational invariance holds 
to a larger extent and the behavior of the finite system is more symmetric 
around the critical point. 
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Figure 1: Lattices PI [6], SH [6] and S[6]. 




Figure 2: The rectangular hatched area as attributed to the site x is defined 
by the dual sites (indicated by crosses) and is proportional to the shaded 
area A^^fX that corresponds to the link from x to x + /i. 
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Figure 3: Specific heat peak values vs. In L for lattices PI , SH , TO and S. 
Here and in all the other figures the error bars are smaller than the symbols. 
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Figure 4: (a) Position of the partition function zero closest to the real axis 
in the complex (3 plane for different lattices. The numbers indicate the base 
length N. The real part of the zeros are closer to their thermodynamic value 
for the sphere-like lattices than for the torus. In (b) we plot the results for 
the spherical surface lattices S. 
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Figure 5: Plot of ln(Im zq) vs. InL; the fit represents the leading FSS 
behavior. 
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Figure 6: Specific heat per unit volume vs. (3 for TO [37] and SH [16]; tliese 
lattices of similar volume are compared to the Onsager solution. One finds 
that the curve for the SH lattice is more symmetric around /3c than for the 
TO. 
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Figure 7: Results for the pseudocritical coupling vs. L for (a) the specific 
heat cy and Vclb, (b) the real part of the closest partition function zero 
and Ui, for pillow, dual cube and torus type lattices. The curves are fits 
according to (|T9|l for the pillow and cubic surface lattices as described in the 
text. 
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Figure 8: Peak positions of the cumulants on spherical surface lattices S. 
One can clearly see, that f3c differs from the critical temperature on the torus 
(^c,TO ~ 0.44068. The curves are from fits to ()19l ) The error bars are smaller 
than the symbols. 



24 



